Stability of Biaxial Nematic Phase in Model Bent-Core Systems 

Piotr Grzybowskj*] and Lech LongaQ 

Marian Smoluchowski Institute of Physics, 
Department of Statistical Physics and Mark Kac Center for Complex Systems Research, 
Jagellonian University, Reymonta 4, Krakow, Poland 
(Dated: November 21, 2009) 

Abstract 

We study a class of models for bent-core molecules using low density version of Local Density 
Functional Theory. Arms of the molecules are modeled using two- and three Gay-Berne (GB) 
interacting units of uniaxial and biaxial symmetry. Dipole-dipole interactions are taken into ac- 
count by placing a dipole moment along the Ci symmetry axis of the molecule. The main aim 
of the study is to identify molecular factors that can help stabilizing the biaxial nematic phase. 
The phase diagrams involving isotropic (I), uniaxial (Njj) and biaxial (Nb) nematic phases are 
determined at given density and dipole strength as function of bent angle. For molecules composed 
of two uniaxial arms a direct / — Nb phase transition is found at a single Landau point, which 
moves towards lower bent angles with increasing dipole magnitude. For the three-segment model 
strengthening of the dipole-dipole interaction results in appearance of a line of Landau points. 
There exists an optimal dipole strength for which this line covers the maximal range of opening 
angles. Interestingly, the inclusion of biaxial GB ellipsoids as building blocks reveals the direct 
I — Nb transitions line even in a non-polar, two-arms model. The line is shifted towards higher 
opening angles as compared to the uniaxial case. 
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Thermotropic biaxial nematic phase, embodying the "holy grail" pj of the liquid crystals 
physics, recently has engendered much scientific interest. In 2004 a long-awaited discovery 
was announced {2, 3], 4\ and the thermotropic biaxial nematics were claimed to be found at 
last, 34 years after its first theoretical prediction by Freiser in 1970 [5|. The discovery was 
made in systems of the bent-core mesogens, also often called banana- or boomerang-shaped. 
They possessed a rigid bent-core with apex angle of ~140°, although recently the biaxial 
phase has also been reported for bent-core systems with apex angle of 90° (|. It should be 
noted that first experimental speculations of possible biaxiality in bent-core systems were 



made earlier 



7]- 



From theoretical point of view it had been known that one can expect a biaxial phase 
in bent-core systems. Bifurcation analysis for hard boomerang model, with molecules com- 
posed of two hard sphero cylinders joined at their ends, was carried out by Teixeira, Masters 
and Mulder [8J predicting phase diagrams with a direct I — N B phase transition at an iso- 
lated Landau point for apex angle of 107°. Subsequent mean-field analysis of interacting 
quadrupoles by Luckhurst 9j have predicted that only a bond angle within a few degrees of 
the tetrahedral angle (109.47°), where Landau point was observed, would allow the biaxial 
nematic to appear above the freezing point of a real uniaxial nematic. That raised a new 
question, since the experiments reported the angle to be near 140°. Among others that is- 
sue has been addressed in various simulation approaches, including Metropolis Monte Carlo 
study of a Lebwohl-Lasher lattice model [lo| . Introduced asymmetry of molecule arms in 
the aforementioned paper shifted Landau point towards lower angles. The same model was 



11 



121 ]. and recently of dipole-dipole 



used to investigate the influence of molecule flexibility 
interactions [13j on phase diagram. 

As it has appeared quite recently, the biaxial phase proves to be more challenging ex- 
perimentally than earlier thought 14J. It seems therefore important to look into molecular 



models that can help in a proper understanding the biaxiality issue. In the present paper 
we construct a class of models of bent-core molecules with the aid of including a broad 
range of molecular factors that can appear relevant in stabilizing biaxial pha se. We build up 
the molecule from parts interacting through Gay-Berne (GB) potential [15| and go beyond 
lattice models to account, at least partly, for excluded volume effects. We will keep the 



arm's length close to experimental values. The simulations for similar molecules^ 



3ut with 



relatively short arms, have not shown any trace of the biaxial nematic phase 16 
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FIG. 1. The construction of pair potential out of two and three GB interacting parts firmly 
attached to molecule. Each part of one molecule interacts with every part from other molecule via 
Vgb- Dipole-Dipole interaction, Vdd, is also added. 

Apart from isotropic and uniaxial nematic phase, smectic phases were observed, which is in 

n 



agreement with predictions o 
The atomistic simulations 



hard sphero cylinder dimer model presented in [20] . 
2l| of the molecular systems studied experimentally 



have 



confirmed the existence of a weak biaxiality in that system. Similarly, for a system of GB 
bent-core molecules with long arms the biaxial nematic phase was observed on cooling [22] . 
Our attempt in the present work is to employ the GB interaction 15j, |23j to model inter- 



acting parts of the molecules. Then, the models are studied using low density approximation 
to the Local Density Functional theory 24, 25]. By means of bifurcation analysis 0, H], 



we determine bifurcation phase diagrams for bent-cores constructed out of two and three 
GB parts, as shown in Fig. HJ Finally, we introduce the dipole-dipole interactions with the 
dipoles taken parallel to molecular C 2 symmetry axis and investigated the influence of the 
dipole strength on stability of Nr. 

We employ reduced units by setting cr = 1 and eo = 1 for the GB potential parameters 
flo| . The reduced distance r* and the reduced temperature t are then given by r* = r / a 
and t = ksT/eo, respectively. We also set u — 1, \i — 2 15j and choose the ratio of length 
to breath of 5:1 for the uniaxial arms. 

The biaxial arms of bent-core molecules are modeled with the help of soft GB ellip- 
soids as proposed by Berardi-Fava-Zannoni 23j, where for the axes we take (a x ,a y ,a z ) = 
(1.2, 0.514, 3.4) and for the potential parameters (e x , e y , e z ) = (1.0, 1.4, 0.2); in addition \i = 1 
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FIG. 2. Exemplary equipotential surfaces for model bent-core molecules composed of 2 ( bottom) 
and 3 (top) uniaxial parts for arm's elongation of 5 : 1 and for opening angle 7 = 126°. The case 
with the dipole-dipole interaction included (p = 2.0) is shown to the right. Surfaces are shown for 
the total pair potential equal to and —0.2. 



and v = 3 (for definitions of the parameters please refer to 23j). The ellipsoids are oriented 
so that the shortest axis is perpendicular to the molecular symmetry plane containing molec- 
ular parts and the longest one is lying in that plane. As mentioned in 23j those parameters 
and orientation of the ellipsoids make the attractive forces strongest in the face-to-face con- 
figuration (that is in the direction of shorter axis). Clearly, that model exhibits molecular 
biaxiality in the limit of 7 = 180° and 7 = 0.0; there should also exist an angle separating 
rod-like and disk-like molecular shapes. 

The resulting GB potential Vqb for a pair of bent-core molecules is then a sum of four 
terms for the two-part case and of nine terms for the three-part case. The dipole-dipole part 
of interaction is of standard form: 

r r ( \ Mi ■ - 3 {hi ■ r) (fi 2 ■ r) 

where /x, = fi£ii,i = 1,2; fa is the unit vector and /x the magnitude of /Xj. The total pair 
interaction, V, is the sum of Vqb an d Vbd- In Fig. [2] shown are the exemplary equipotential 
surfaces. Figures show the profiles for four different relative molecular orientations, including 
two- and three GB part bent-core molecules of uniaxial arms. 
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TABLE I. Dipole-dipole contribution to total pair potential energy in ground state. 



H \Vdd/{V G b + V dd )\ 

2.8 0.50 

2.2 0.40 

1.6 0.25 

1.5 0.22 

1.2 0.15 



Relative importance of Vqb and Vdd parts is measured by by the ratio 



v Dl 



The 



Vgb+Vdd 

ground state values of this parameter for different values of p. are given in Table [H We study 
the equilibrium properties of the systems by minimizing the grand potential with respect to 

The necessary condition is given by 



the one-particle distribution function p |24j . 



25 



26, 2 



a self-consistent integral equation for stationary distribution ps(q)- It reads [27 ] 



Ps(q) 



Z s l exp{Ci(g, \p s ])}, 



(1) 



where Z$ = f exp{Ci(g, [ps])}dq/ (N) is the normalization constant, Ci(q, [ps]) is the one- 
particle direct correlation function, (N) is the average number of particles in the system and 
q represents collectively all degrees of freedom of a single molecule. 

In the present study we (a) approximate C\ using second order virial expansion, some- 
times also referred to as low density approximation and (b) restrict analysis to nematics, 
which amounts in setting ps(q) = pP(fi), where p = < y > &o is the dimensionless density and 
f2 are the Euler angles parameterizing orientation of a molecule-fixed frame with respect to 
the laboratory-fixed frame. After these limitations the Eq. ([I]) becomes: 



p(n x ) = z- 1 



cxp 



p / C 2 (f2 1 f2 2 )P(fi 2 )df2 2 



(2) 



with 

c 2 (fiiS7 2 ) 
and normalization constant: 



exp 



—v(n^n 2 ,r* 12 ) 



1 > d 3 r* 2 , 



exp 



p / C2 (o 1 n 2 )p(rj 1 )do 1 



dfi 2 . 
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Here dO = dad (cos (/?)) oVy stands for integration over Euler angles, d 3 r* = r* 2 dr*d (cos (9)) dip 
and 0^ x 2 is the relative orientation of the molecules. 

The usual expansions of -P(O) in the base of D 2h symmetry adapted A functions [26] 

p w = E ^^Ma&w (3) 

L,m,n 



allows to introduce the order parameters Amjn — J dOP(O) A^i(n) for nematics. The 
summation runs over the allowed values of {L,m,n}, where L is a non- negative integer. 
Generally, if L is even, then < m < L and < n < L. If L is odd, then 2 < m < L and 
2 < n < L. If, in addition, we expand c 2 : 

^(OxOs) = J2 C mnA£UOr 1 2 ) 

L,m,n 

where 

2L + 1 



do^o 2 c 2 (0!0 2 ) A^j n (n^n 2 ) 

and where c mn = c nm due to particle interchange symmetry, then Eq. (2) becomes reduced to 
the set (in general infinite) of nonlinear equations for the order parameters. Using bifurcation 
analysis we now seek for a subset of nonzero order parameters, describing low-symmetry 
phase, that branch off from the background high-symmetry phase. Generally, in the isotropic 
phase all order parameters vanish. The uniaxial phase is characterized by nonzero order 
parameters indexed by m = 0. Finally, in the biaxial nematic phase all the order parameters 
become nonzero. 

The bifurcation points so determined are either spinodal points for the first-order phase 
transitions or critical points for the continuous transitions. Hence, for continuous and weakly 
first-order phase transitions, as holds of isotropic and nematic phases, we arrive at quite ac- 
curate estimates of the phase diagrams in temperature-density plane. Following the analysis 
as described in [23] two different bifurcation formulas can derived. The first one is the 
equation for the bifurcation from the isotropic phase. It reads 

P = / • ( 4 ) 



coo + c 22 - y 4c5 2 (c 00 - c 22 ) 
For the uniaxial to biaxial bifurcation a more complex formula is found 

2c o (adc 02 - d) + c 2 (2a + a6c 2 2) - d 2 c o + bc 2 2 



p = 35 



(a 2 + 2bd) (cg 2 + C00C22) 
a/ a 2 c 00 C22 + bd (c 00 C22 - 2c 2 , 



^02) 



(a 2 + 2bd) {c 2 02 + c 00 c 22 ) 



(5) 
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FIG. 3. Diagrams for two-parts banana, for two densities (p) and three values of dipole moment 
(//). On each plot two branches of bifurcation from uniaxial to biaxial phase meet the upper line 
of bifurcations from isotropic phase. 

where 



a = 20 A^ + VTE A ( qI 



b = 14 + 20A< 2 > + Aqq + V35A^ 



(4) 



d=7- 10Ag> + 3A^ , 



.(4) 



and where Amjn are determined in the uniaxial nematic phase. We would like to point 
out that since the formulas for Am,'„s depend on p, the equation (jSJ) for given t becomes a 
self-consistent equation for density. 

Please note that each CfYin IS Si SIX dimensional integral. We performed numerical in- 
tegration to obtain temperature dependence for each of the coefficients for given set of 
the molecular parameters. The integration procedure was implemented in C. We have in- 
corporated both Monte Carlo and adaptive multi-point Gauss quadratures method in the 
integration procedure. We performed Monte Carlo (MC) integration over orientations us- 
ing quaternion parametrization of rotations and then calculated the integral over length of 
intermolecular vector ri2 using adaptive Gauss quadratures. Approximately 2 million of 
MC cycles were used to calculate the integral; the relative error was estimated to be less 
than 1%. Now we present the bifurcation diagrams, which follow from the solutions of 
the bifurcation equations (jlj) and (JSJ for the uniaxial Vqb- In Fig- [3] and H] we have plotted 
t divided by t 90 , the temperature of bifurcation from isotropic phase for 7 = 90° in the 
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FIG. 4. Diagrams for three-parts banana, for two densities, and five dipole moments. Each plot 
shows the bifurcation from isotropic phase and two lower branches of uniaxial-to-biaxial phase 
transition. 



TABLE II. Bifurcation temperatures from isotropic phase for 7 = 90°. 
P 

bent-core molecule composed of two parts: 
0.04 
0.06 

bent-core molecule composed of three parts: 
0.026 
0.04 



*90 

0.84 
0.96 

1.34 
1.71 



non-polar case. Numerical values of tgo are gathered in Table [TTl Note that the standard 
phase sequence jsj is recovered that involves the isotropic phase, the rod-like and disc-like 
uniaxial nematic phase and the biaxial nematic phase. For the case of two-part molecule 
without the dipole moment Landau point is found to be near 7 = 107°, in agreement with 
the hard-boomerang model The diagrams include two dipole strengths /x = 1.2 and 
\i = 1.5 for which dipole-dipole interaction is about 15% and 22%, respectively, of the total 
potential in ground state. The bicritical point is shifted towards lower angles with increasing 
/i, which contrasts with the results of Monte Carlo simulation for the Lebwohl-Lasher lat- 



tice model 



131 ]. where introduction of the dipoles resulted in a line of direct isotropic-biaxial 
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TABLE III. Landau point versus dipole magnitude /x 



/x 

bent-core molecule composed of two parts: 
0.0 
1.2 
1.5 

bent-core molecule composed of three parts: 

0.0 
1.2 
1.6 
2.0 
2.8 

bent-core molecule composed of two biaxial ellipsoids 

0.0 



Landau point 

107° 
104° 
103° 

p=0.026 /9=0.04 

89° 89° 

86° 86° 
74° - 86° 74° - 86° 
63° - 86° 63° - 80° 
83° - 97° 82° - 92° 

121° - 128° 



transitions. 

For the non-polar three-part molecule, the Landau point is found to be at 7 = 89° and 
is shifting to lower angles with increasing dipole magnitude (Table HIT]) up to a point where 
the dipole-dipole interactions make up 20% (/x = 1.4) of the total potential in the ground 
state. Above that value the bicritical point changes into a line of Landau points that widens 
with increasing /x; for /x = 1.6 it covers the range of 12° and for /x = 2.0 it extends for more 
than 20°. 

The low 7 boundary practically does not change (for lower density) and is equal to 86° 
for the dipole strength /x < 2.1. Then the bicritical region begins to shrink and is shifted 
towards higher angles. The highest dipole studied was the one of /x = 2.8 for which Vdd 
approaches 50% of the total potential energy. As can be seen from Fig. [3] the bicritical line 
in that case is still getting shorter and moves towards higher bond angles. Fig. [6] shows the 
evolution of Landau region as function of the dipole magnitude, /x. 

The diagrams are presented for two densities such that the corresponding packing fraction 
is of the order of 0.3—0.4. As can be seen from Table UTTl some differences appear with varying 
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FIG. 5. Bifurcation diagram for non-polar bent-core molecules modeled by two biaxial GB ellip- 
soids \2c 



density for the strongest dipoles (fi > 2.0). Namely the line of the direct isotropic-biaxial 
transitions shrinks for higher density. 

Finally we take into account the model where the arms of the molecules are biaxial. We 
are going to address the issue of observed disagreement between the angles for which the 
Landau point appears experimentally (7 = 140°) and theoretically. Results presented below 
replace uniaxial GB arms with their biaxial version developed by Fava, Berardi and Zannoni 



231 ] . Interestingly, the Landau point in the biaxial model, Fig. [5j is replaced by a line of 
bicritical points even for the non-polar molecule. That line starts near 121° and ends for 
7 = 128°. The region becomes reduced to a single point with decreasing arm's biaxiality. 

Summarizing, we have presented a bifurcation study for a class of models with character- 
istic features typical of the bent-core molecule. Using Density Functional Theory we have 
retrieved the diagrams in low density approximation. Analysis included two and three-part 
bend-cores with arms modeled by GB interacting ellipsoids of uniaxial and biaxial sym- 
metry. The dipole-dipole interaction was added and the dipole strength influence studied. 
Non-polar uniaxial model revealed a single Landau point, in agreement with results for hard 
molecules js]. The deviation from uniaxial symmetry of the arms resulted in transformation 
of the single bicritical point into a line of direct isotropic-biaxial transitions. The inclusion 
of the dipole-dipole interactions resulted in shifting of the Landau point towards lower bond 
angles with increasing dipole magnitude. For the case of two-arm molecule stronger dipoles 
easier destabilized uniaxial nematic phase. For the three-part banana a line of the bicriti- 



10 



95° 
90° 
85° 

c 
o 

'§? 80° 

c 

I 75° 
(D 

- 1 70° 
65° 
60° 

0% 10% 20% 30% 40% 50% 

Dipole-dipole contribution in total potential 

FIG. 6. The range of direct isotropic-biaxial transitions in bend angle, in function of dipole-dipole 
contribution in total potential. 

cal points has emerged. The results suggest that there exists an optimal dipole range that 
makes the appearance of the biaxial phase most probable. 
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